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ABSTRACT 

A loss of magnetic field confinement is believed to be the cause of coronal mass 
ejections (CMEs), a major form of solar activity in the corona. The mechanisms 
for magnetic energy storage are crucial in understanding how a field may possess 
enough free energy to overcome the Aly limit and open up. Previously, we have 
pointed out that the accumulation of magnetic helicity in the corona plays a 
significant role in storing magnetic energy. In this paper, we investigate another 
hydromagnetic consequence of magnetic-helicity accumulation. We propose a 
conjecture that there is an upper bound on the total magnetic helicity that a 
force-free field can contain. This is directly related to the hydromagnetic property 
that force-free fields in unbounded space have to be self-confining. Although a 
mathematical proof of this conjecture for any field configuration is formidable, 
its plausibility can be demonstrated with the properties of several families of 
power-law, axisymmetric force-free fields. We put forth mathematical evidence, 
as well as numerical, indicating that an upper bound on the magnetic helicity 
may exist for such fields. Thus, the accumulation of magnetic helicity in excess of 
this upper bound would initiate a non-equilibrium situation, resulting in a CME 
expulsion as a natural product of coronal evolution. 

Subject headings: MHD — Sun: magnetic fields — Sun: corona — Sun: coronal 
mass ejections (CMEs) 
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1. Introduction 

Coronal mass ejections (CMEs) are a major form of solar activity in the corona. A CME 
takes 10^^"^^ g of plasma from the low corona into the solar wind to disturb the near-Earth 
space if the CME propagation is directed toward us. UnUke other solar activities such as 
flares and filament eruptions, CMEs were only first observed when white light coronagraphs 
were put into space in the 1970s (MacQueen et al. 1974, Howard 1976). Since then, this 
astrophysical process has caught much attention in the space physics community because of 
the adverse consequences of CMEs for orbiting satellites. 

CMEs are identified as large-scale bright features in white hght that move outward 
through the corona at speeds from 10 to over 2000 km/s (Hundhausen 1999). These are 
impressive dynamical phenomena that we happen to be able to observe with a certain degree 
of completeness in both white light (Dere et al. 1999, St. Cyr et al. 2000, Burkepile et al. 
2004) as well as in X-ray and EUV bands (Webb & Hundhausen 1987, Hudson & Cliver 
2001). In a low-/? plasma system on the large scales of the corona, we generally believe 
that MHD processes control the dynamical evolution of the magnetic field and the coronal 
plasma. In considering CMEs as a MHD phenomena, previous studies have mainly focused 
on two physical issues: magnetic energy storage and CME initiation. 

The energy storage problem is concerned with how much free magnetic energy can be 
stored in the corona in order for the field to spontaneously open up to infinity during a 
CME eruption. This issue was first brought up in Aly (1984), and later on in Aly (1991) and 
Sturrock (1991), where Aly conjectured that a fully- anchored force-free field may not be able 
to store enough free magnetic energy to drive a CME kinematically, in addition to opening 
up the field under conditions of high electrical conductivity. This conjecture states that an 
anchored force- free field cannot store more energy than that of any fully-open field sharing 
the same boundary fiux distribution at the base of the atmosphere, where those fully-open 
fields have energies bounded from below by a minimum usually referred to as the Aly limit. 
Many studies have pointed out that the constraint of the Aly Conjecture can be bypassed 
by relieving one of the conditions assumed in the Conjecture. For example, by relieving 
the force-free condition an equilibrium field with plasma weight can easily store enough free 
magnetic energy to drive the CME (Low & Smith 1993, Wolfson & Dlamini 1997, Fong et 
al. 2002, Low et al. 2003, Zhang & Low 2004, Flyer et al. 2005). Even under the force-free 
condition, recent studies have showed that the Aly hmit can be exceeded if the field contains 
a detached magnetic flux rope in the corona (Hu et al. 2003, Wolfson 2003, Flyer et al. 
2004). The Aly Conjecture can also be bypassed through the topology of multipolar fields. 
For multipolar fields, a CME- like expulsion only needs to open up one but not all the bipolar 
fields (Antiochos et al. 1999). In such a situation, the Aly limit associated with opening up 
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the entire field is not the relevant energy threshold; see Zhang & Low (2001) for illustrative 
examples. 

The CME initiation problem is concerned with the sporadic behavior and locations of 
CME occurrences. CMEs are occurring at an average rate of one to three events per day, with 
the frequency of occurrence increasing as solar activity intensifies. The catastrophe model 
(Forbes & Isenberg 1991, Isenberg, Forbes & Demoulin 1993) interprets the CME initiation 
as a catastrophic loss of equilibrium in a series of quasi-static evolution of magnetic equilibria 
corresponding to mathematical variations of some model parameters. The fiux-emergence 
trigger model (Chen & Shibata 2000) proposes that the newly-emerged flux near a preexisting 
flux rope is the trigger of the expulsion. 

It might be useful to point out that when we address the magnetic energy storage 
problem or the CME initiation problem, our underlying thinking tends to take these as two 
independent issues. However, as we have pointed out in a recent review (Zhang & Low 2005) 
these two issues can be related by the accumulation of magnetic hehcity in the corona as a 
fundamental MHD process. The accumulation of magnetic helicity in the corona leads to 
a magnetic energy build-up as a natural product of coronal evolution. Eruptions, as a loss 
of fleld conflnement, occur when the stored magnetic energy exceeds the Aly limit. In this 
paper, we take another approach to address the role of magnetic helicity accumulation in 
the corona. We propose a hydromagnetic conjecture that there is an upper bound on the 
total magnetic helicity that a force-free field can contain, owing to the fact that a force- 
free field in an open atmosphere has to self-confine. This conjecture poses a fundamental 
but mathematically challenging problem in magnetohydrodynamics. We focus our attention 
in this paper on demonstrating the physical plausibility of the conjecture by using several 
families of numerical solutions of nonlinear force-free fields to make our physical points. 

The paper is organized as follows: In §2, the physical considerations on the roles of mag- 
netic helicity accumulation are addressed, with a particular focus on their connections to the 
magnetic energy storage and CME initiation. In §3, power-law axisymmetric force-free flelds 
are used to indicate the possible existence of an upper bound on the total magnetic helicity, 
using both analytical inequalities and numerical solutions. Conclusions and discussions are 
presented in §4. 
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2. The Physical Picture and the Conjecture 

2.1. Magnetic helicity and its approximate conservation in the corona 

Magnetic helicity is defined as H{V) = JyU ■ 'BdV, where u; is the vector potential of 
the vector magnetic field B in a volume V. It is a physical quantity that quantitatively 
measures the topological complexity of a magnetic field, such as the degree of linkage and/or 
twistedness in the field, provided V wholly contains B. We henceforth refer to F as a a 
magnetic volume. With this physical quantity, the frozen-in condition in a perfectly con- 
ducting plasma can then be stated in terms of the absolute conservation of magnetic helicity 
in all magnetic sub- volumes as well as of the total magnetic helicity for the whole magnetic 
volume (MofFatt 1985, see Zhang & Low 2005 for more descriptions). 

In a situation of a domain V whose magnetic field B threads across its boundary, i.e., 
the normal magnetic field Bn does not vanish across the domain boundary dV, H{V) is then 
gauge dependent and is therefore not physical. Magnetic twist and linkage in V then must 
be described by the relative magnetic hehcity denoted by Hji{V) introduced by Berger Sz 
Field (1984). We will define Hr{V) later with the development of the paper. In this paper, 
we take magnetic helicity to mean either HiV) or Hjiiy), depending on whichever applies 
to the given situation. 

The frozen-in condition, or the absolute conservation of magnetic helicity in magnetic 
volumes, must break down in a realistic plasma where the electrical conductivity is very 
high but not infinite. Electric current sheets can form, and magnetic reconnection will 
subsequently set in, to dissipate the thin current sheets and change field topology (Parker 
1994). However, by a general scaling analysis of resistive processes, Berger (1984) pointed 
out that the dissipation rate of the total relative magnetic helicity in the corona is orders of 
magnitude smaller than the release rate of magnetic energy. So the total magnetic helicity in 
the corona is still approximately conserved despite the occurrences of magnetic reconnection. 

The conservation of total magnetic helicity in the corona has a profound influence on 
the coronal dynamics. Magnetic helicity is continuously being transported into the corona 
(Kurokawa 1987, Lites et al. 1995, Leka et al. 1996, Zhang 2001, Pevtsov, Maleev & Long- 
cope 2003) with a preferred sign for the two solar hemispheres, positive and negative in the 
southern and northern hemisphere, respectively (Rust 1994, Pevtsov, Canfield & Metcalf 
1995, Bao & Zhang 1998). Moreover, the same opposite signs of heficity are statistically pre- 
ferred in the two respective hemispheres, negative and positive, respectively, in the northern 
and southern hemispheres, independent of the solar cycle. Thus, under the approximate 
conservation of total magnetic helicity, the hemispherical preferred signs of helicity implies 
an accumulation of magnetic helicity in the two hemispheres. It is this accumulated helicity 
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that, we propose, characterizes the build up of magnetic-energy and when accumulated to a 
point a magnetic field must open up as an eruption. 

2.2. Helicity-energy relationship for fields in finite domains 

For a force-free field confined in a finite domain with rigid walls, its total magnetic 
energy can exceed any value by endowing the field with a sufficiently twisted topology. This 
can be easily seen from the Virial Theorem of Chandrasekhar (1961), which states that for 
a force- free field confined in a sphere of radius r = 1, its total magnetic energy, given in 
spherical coordinates, is 




For a fixed boundary flux &t r = 1, the more highly twisted the force- free fleld is, the 
greater are the boundary values of the tangential components Be and B^p at r = 1, leading 
to a higher total magnetic energy of the field. Note that force-free solutions can exist in such 
a domain even if 5^ = on r = 1. 

More importantly, if the field is given a total non-zero magnetic helicity Hq, Woltjer 
(1958) Theorem sets a ground-level energy above that of the potential field. The Woltjer 
Theorem produces many constant a force-free fields, each being a local minimum-energy state 
compared to neighboring flelds with the same total magnetic helicity Hq. The absolutely 
lowest minimum-energy state, taken from among all those constant a force-free fields, is 
called the Woltjer state, characterized by a specific constant ao. This constant ckq of the 
Woltjer state, together with the given total magnetic hehcity Hq, puts a bound from below 
for the magnetic energy of the field in a finite domain. That is, the magnetic energy of any 
field B possessing the same total magnetic helicity of Hq satisfies the bound: 

E > I uqHq I . (2) 

Applying this theorem to a coronal magnetic field underlying a locally confined coronal 
structure (Zhang & Low 2003), we can see that the coronal field cannot relax to a potential 
state as long as the total magnetic helicity is not zero. Magnetic reconnection may set in to 
release the excess magnetic energy that the coronal field may happen to have, but the field 
cannot relax to a state that possesses less magnetic energy than that of the Woltjer state. 
With more and more magnetic helicity brought with a preferred sign into the corona, this 
ground-level magnetic energy increases with the increasingly accumulated helicity. This is 
the principal consideration described in Zhang & Low (2005) of the role of the accumulated 
magnetic helicity in storing magnetic energy for CMEs. 



6 



2.3. 



Confinement issue for fields in unbounded domain 



For a forcc-frce field filling an unbounded atmosphere, it has an energy property that is 
distinct from the one confined in a finite domain with rigid walls. Now, the Virial Theorem 
of Chandrasekhar (1961) takes a different form and sets a special condition for those fields 
in force-free equilibrium. For a force-free field outside a sphere of radius r = 1 in spherical 
coordinates, its total magnetic energy now takes the form of 



where we assume that |B| spatially decline to zero at r — > oo so that the total magnetic energy 
Eq is finite. Here and elsewhere we assume this is the case whenever we treat integration 
over unbounded domain without explicitly stating this qualification. Since Eq should always 
be positive, the above equation says that no force-free state is possible in r > 1 domain 
unless Br on r — 1. With Br 7^ 0, equation (3) sets an absolute upper bound on the 
total magnetic energy. 



which is uniquely determined by the flux distribution on r = 1, independent of the twist of 



The essence of this hydromagnetic result is that there is an upper bound on the total 
magnetic energy that a force-free field can attain and this upper bound is determined by its 
boundary fiux at r = 1. This is equivalent to saying that a force-free field in an unbounded 
atmosphere requires certain means of confinement to remain in equilibrium. The field with 
with an arbitrarily large, prescribed degree of twist may not be able to find a force-free 
state to meet the bound on its total magnetic energy. This is indicated by the the following 
analysis. 

A field completely detached from r = 1, as drawn in Figure lb where Br = 0, can never 
be in equilibrium in the absence of other body forces. The only way to hold a flux rope, and 
hence attain a certain amount of magnetic helicity in the low corona, is to add a 7^ 
fleld as in Figure la to function as an anchoring agent. The resultant conflguration will be 
the one sketched in Figure Ic, a fleld containing a certain amount of magnetic helicity and 
yet flnding equilibrium in the r > 1 domain. 

Now imagine we flx the Br distribution on r = 1 in Figure Ic but arbitrarily increase 
the total azimuthal flux in the flux rope. This parametric change will take the fleld closer 
and closer to the conflguration in Figure lb. We intuitively expect that there would be 
a parameter representing the fleld topology crossing over where no force-free equilibrium 




(3) 




(4) 



the fleld. 
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can exist. We suggest that this critical parameter corresponds to a point of maximum 
accumulation of magnetic helicity in the atmosphere. 

This conjecture intuitively comes from the understanding that the total magnetic helicity 
cannot be destroyed in the corona at observable time scales. So the corona may thus get a 
chance to accumulate magnetic helicity of a preferred sign to cross over a threshold bounding 
the allowable total helicity in a force-free field. A rigorous mathematical bound on the total 
helicity of a force- free field in r > 1, in full generality, is needed to prove the conjecture, 
in parallel to equation (3) for the bound on the total energy of a force-free field. This is a 
difficult task in mathematical physics. As a first step towards such a proof, we use families 
of power-law axisymmetric force-free fields of Flyer et al. (2004) to investigate the nature of 
this conjecture. We focus the rest of our paper on such a study in order to provide a basic 
understanding of the physical issues. 



3. Helicity Bounds for Power- law Axisymmetric Force- free Fields 
3.1. A Formula of Total Relative Magnetic Helicity of Axisymmetric Fields 



We first derive a formula to calculate the total relative magnetic helicity in axisymmetric 
magnetic fields. This simple formula is limited to axisymmetric fields but is otherwise quite 
general, independent of whether the axisymmetric field is in equilibrium or not. We therefore 
derive the formula from first principles. 

Under axisymmetry, we can always write the solenoidal magnetic field B in r > 1 in the 
form of 



B = 



r sin 6 



IdA _ dA ■ 
r 09 ' Or '' 



(5) 



Here, the fiux function A defines the poloidal magnetic field and the function Q defines the 
toroidal field. 

Alternatively, we can also use a vector potential uj to define a solenoidal magnetic field, 
where the magnetic field B can be derived from the vector potential uu by the following 
relationship: 

B = V X a; . (6) 
Then for an axisymmetric magnetic field, uj takes a form of 



1 

rsinO 



'IdQ' _dQ' 
r 86 ^ dr 



(7) 
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where Q' and A! are two physical quantities that are related to the functions A and Q by 
with 11 = cos 9, and 

A' ^ A + constant . (9) 

These relationships between Q', A', Q and A can be understood by the following analysis. 
From the definition of the vector potential and Stokes Theorem, the magnetic fiux across any 
surface S with a boundary C is F = J^B ■ dS = J^,uj-d\ where d\ is a path element along C. 
Therefore, the vector potential a; carries information about the magnetic flux across surfaces. 
Application of this relation with uj given by equation (7) for an axisymmetric field shows 
that the function Q' carries information on the azimuthal fiux, as implied by the relationship 
between Q' and Q in equation (8), and the function A' carries information about the poloidal 
fiux as expressed by equation (9). 

The constant in equation (9) can be used to take care of the physical requirement that 
A' should be zero along the polar axes. Without loss of generality, we set this constant to 
zero, which is equivalent to setting the fiux function A to zero along the polar axes. 

Suppose we are given the field B in the unbounded space r > 1 denoted as V. This field 
has a fiux anchored to the inner boundary r = 1 so that the classical hehcity of Woltjer is 
not gauge invariant and we must use instead the relative hehcity of Berger & Field (1984). 
This is given by the formula derived in Berger (1985): 



Hr^ / ujP-Bl dV+ / uj-B dV , (10) 

JVa JV 

which involves extending the domain V to the spherical region r < 1 denoted as Va where 
is a potential field with vector potential a;^. This potential field continues into the given 
field B across the boundary with continuity of Br- The important point here is that Hr is 
independent of the gauge and is a physical quantity describing the topological complexity of 
the field B in r > 1. We calculate Hr as follows. 

The solenoidal condition requires 

B^r |,=i = B-r (11) 

for the continuity of the normal field across r — 1. Equivalently, in terms of the vector 
potential, we require 

X r |r=i = a; X r |r=i , (12) 
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for the continuity of the tangential components of the vector potential (Berger & Field 1984). 

Equation (12) requires that both A and ^ are continuous across r — 1. The continuity 
of A across r = 1 is just equation (11) and the normal flux of B at r = 1 determines a unique 
potential field in 14- This potential field is generated by a flux function A^, with Q — 0. 
We take Q' — everywhere in T4, which is consistent with equation (8) where Q — 0, as a, 
selected gauge. Hence, we have 



IdAl 



rsinO \r 89 
A ^ 



r — 



or 



rsinO' 



(13) 



which are two orthogonal vectors so that Jy u;^ • ^^dV — and Hji in equation (10) then 
contains only the contribution from r > 1. 



Now we have the total relative magnetic helicity in r > 1 as 

1 



H 



R 



j 

Jr>l 



uj ■ BdV = 27r 



[VQ' -VA + AQjr^drdu , 



(14) 



where we have the geometric boundary conditions that both Q' and A vanish at the polar 
axes and the gradients of these two functions with distance vanish sufficiently fast at infinity. 
These conditions and integrations by parts allow us to rewrite 



1 



-Vg' • \/Ar^drdfx 



drd/i- 



AQ 



'r>l r^C^ - )"^) ' " ""' " ' Jr>l ""1 - 

And then the total relative magnetic helicity is given by: 

-1 



1-/.2 



A 



dr 



(15) 



r=l 



H 



R 



Att I drdfj, 

ir>l 



= 47r 



/ 

Jr>l 



drd/i 



AQ 

1-/^2 

AQ 



271 



dfi 



1-/.2 



A 



dQ[ 
dr 



r=l 



a simple form of an integral to calculate the total relative magnetic helicity in r > 1. 
Note this formula can also be rewritten as 



(16) 



H 



R 



An 



A X B^rdrdO 



(17) 



and it gives a physical meaning of this simple form of helicity formula: The total relative 
magnetic helicity in axisymmetric fields is just a simple convolution of the local wrapping of 
the azimuthal fiux element B^rdrd9 with the poloidal flux A. 
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3.2. Power- law axisymmetric force- free fields 

We adopt the families of numerical force-free fields from Flyer et al. (2004) where Q as 
a strict function of A has the form of 

Q^i^A) = ■ (18) 

Here n is an odd constant index required to be no less than 5 in order for the field to possess 
finite magnetic energy in r > 1 and 7 is a free parameter which we choose to be positive 
without loss of generality. We refer this family of force-free fields power-law axisymmetric 
force-free fields in this paper. 

Taking this form of Q, the force-free condition, that is, 

(V X B) X B = , (19) 
reduces to the following governing equation for flux function A: 

This governing equation was solved numerically as a boundary value problem within domain 
r > 1 in Flyer et al. (2004), subject to the prescribed boundary flux distribution 

A\r=i^sm^e . (21) 

We refer interested readers to that paper for various properties of this family of force-free 
flelds for the cases of n = 5, 7, 9. 

Note that although 7 > is a free parameter, it actually has a least upper bound 70 
above which no solution of the boundary value problem can exist. As shown in Flyer et al. 
(2004), with the boundary condition (21), for each flxed n, 

2(n+l)(2n + l)n 

3(2n)!! ^ ' 



3.3. A bound on the total magnetic helicity for a fixed n 

The Virial Theorem expresses the total magnetic energy of a force-free magnetic field 
B in r > 1 in terms of a surface integral given by equation (3). Aly (1988) presented a 
"generalized scalar virial equation" that expresses the total magnetic energy in other forms. 
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Writing S| — 3^+ B'^ and introducing an arbitrary function /(r) with derivative /'(r), Aly 
derived 

/ [f{B'^+B^)+rf'iB^-B'^)]dV = f{r = l) [ {Bl - Bl)dS . (23) 

Jr>l Jr=l 

For our purpose, we set f — r to obtain 

/ [r{B', + B^) + r{B^ - B',)] dV ^ I {B^ - B^)dS . (24) 

Jr>l Jr=l 

This leads to an alternative expression of total magnetic energy of force-free fields as 

Eo ^ ^ [ r {Bl + Bl) dV . (25) 
Prom this form of £'o, we then have 





> 


1 

An 


Eq 


> 


1 

An 



(26) 



r>l 



and 



rB^dV (27) 

r>l 



El > f ^-^dV f '^dV . (28) 
Jr>i 47r y^>i An 

Consider Cauchy-Schwartz inequality 

J^dx j^g^dx>(^jjgdx^ , (29) 

which relates two functions / and g in space X, where dx is the spatial integration element. 
Applying Cauchy-Schwartz inequality to inequality (28), we get 

El > (j^ f rBoB^dv) . (30) 



The right hand side of this inequality contains an integral of the product of the poloidal field 
Bff with the toroidal field B^, suggestive of the formula (16) of the relative magnetic helicity. 

By equation (16) the total relative magnetic helicity in r > 1 of a power-law axisym- 
metric force-free field is 



27 r dn 
nTl 1^2 



00 

n+3 



Hr = 47r;/^^ / / A—dr . (31) 
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Evaluate the right hand side of inequahty (30) we obtain 

1 dfx . 2 



rar-:^\l -A 2 



dr V + 1 



2 J_i 1 - /i^ + 3 V n + 1 
1 da 



27 ■ n+3 , 

-A 2 L=i 



CO 2 

dr 



n + 3 V n + 1 



27 n + 3 



(32) 



with an integration by part. Using the boundary condition at r = 1 to evaluate the first 
term and denoting it as ^^-^/^^Cq, then 



1 da n+3, 

is a constant independent of 7. Inequality (30) now takes the form of 

2 

ri + 3 ^ I I 27 ^ , 



and we have achieved a bound on Hr. 

A simple analysis reduces inequality (34) to the form of 

3 



(33) 



(34) 



-Stt 



n 



-Eo + 



27 
n + 



-Co 



<HR<87r 



n + 3 



En — 



27 

n + \ 



-Co 



(35) 



Since Cq is defined strictly by the boundary condition, inequality (35) shows that the mag- 
nitude of Hr is bounded from above by bounds that are fixed for each given n. 

The above inequality shows that \Hr\ is bounded for each fixed n. However, this does 
not assure us that a single bound Hmax > exists such that \Hr\ < Hmax for all n, because 
the upper and lower bounds go to ±00 as n — > 00. A tighter bound on Hr is then motivated 
to be found and will be presented in §3.5. 



3.4. Numerical evidence 

It is useful for us to use those numerical solutions in Flyer et al. (2004) to examine 
some helicity properties of this family of power-law axisymmetric force-free fields. In Flyer 
et al. (2004), Newton's iteration combined with a pseudo-arc length continuation scheme 
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was used to solve the nonlinear partial differential equation (20) in the unbounded domain 

r > 1 for the boundary condition (21). It was found that iterations converge to solutions only 
for a finite range of 7 for a fixed index n, as indicated by inequality (22). In addition, the 
solution A was found to be a multi- valued function of 7, i.e. for a given value of 7 multiple 
solutions can exist. These multiple functions can be pieced into a single solution curve as a 
function of 7 joined at turning points. The pseudo-arc length continuation scheme was used 
to march along the solution curve for a given n and to pass through the turning points in 
the parameter space where the Jacobian of the Newton's iteration become singular. 

Each point on the obtained sohition curves describes a power-law axisymmetric force- 
free field in r > 1 with a total magnetic energy Eq, a total azimuthal fiux 



and a total relative magnetic helicity Hr given by equation (16). Flyer et al. (2004) char- 
acterized their solutions in terms of £"0 and F^, but not H^. The approximate conservation 
law for naturally suggests the use of H^i to characterize the twisted state of the force-free 
field in r > 1. In Figure 2 we display the variation of Eq and F^ with 7 for the n = 5, 7, 9 
solution curves, taken from Flyer et al. (2004), together with the respective variations of Hji 
for these curves. 

For each of the n = 5,7,9 solution curves, Flyer et al. (2004) found that F^p is bounded 
to be of the order of twice the total poloidal flux given by boundary condition (21). This 
numerical evidence indicates that force-free fields can only self-confine its magnetic pressure 
by its tension force provided its azimuthal fiux is less than a bound determined by its poloidal 
fiux fixed by the boundary condition. These are indicated in Figure 2 by the azimuthal fiux 
Ftp monotonically increases from to its maximum value along a solution curve, noting in 
particular that the maximum value of F^ is in each case less than 2 and not sensitive to n. 
The corresponding relative helicity H^i shows interesting undulations superposed on a general 
trend of monotonic increase along the solution curves, showing a bound on Hf> less than 15 
and not sensitive to n. These numerical solutions provide a basis for the suggestion that 
the total relative magnetic hehcity of the three famihes of power-law axisymmetric force-free 
fields is also bounded by a suitably defined bound that depends only on the boundary fiux 
distribution. 

Figure 3 presents these solutions again, but are plotted for the total magnetic helicity 
(left panels) or for the total magnetic helicity divided by the total magnetic energy (right 
panels) against the azimuthal fiux. These plots show that the total magnetic helicity increases 
nearly monotonically with the increasing azimuthal fiux, which confirms our intuition that 
the confinement issue discussed in terms of the azimuthal fiux is not much different from the 




(36) 
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discussion in terms of magnetic helicity. However, the total magnetic helicity is known to 
be conserved in the corona and therefore accumulated in the corona, whereas the physics is 
less specific whether the azimuthal flux is conserved in the corona or not. 

Another interesting observation from Figure 3 is that, the ratio between the total mag- 
netic helicity and the total magnetic energy shows a better monotonic variation with the 
azimuthal flux than does the magnetic helicity itself. This may reflect a general rela- 
tionship among the three quantities. Notice that inequality (35) suggests that H^i/Eq is 
roughly smaller than 47r(n + 3) which is confirmed by the numerical solutions. However, 
the bound obeyed by the numerical solutions is a lower or tighter upper bound. That is, 
Hr/Eq < Qtt < 47r(n + 3) = (32 — 48)7r, the last term taken from inequality (35). This hints 
that the upper bound we have found by inequality (35) is not a stringent bound and that 
lower upper bounds may exist. 



3.5. Inequality relating total magnetic helicity to azimuthal flux and magnetic 

energy 

The numerical relationships among Eq, E^p and Hr in Figure 3 suggest that a rigorous 
relationship may exist between these three physical quantities, at least, for the power- law 
force-free fields. This rigorous relationship is derived in this subsection by first introducing 



E, = lljliV, (37) 



which is the contribution of the azimuthal field component to the total magnetic energy. 

As an extension of Cauchy-Schwartz inequality. Holder inequality states that, given any 
two functions / and g in space X with spatial integration element dx, we have 

x^J \Jx^''j -Jx^^'^''' ^^^^ 



where k > 1 and is related to k' by 



11 



The inequality (38) reverses direction if A; < 1 while k still relates to k' by equation (39). 
Applying Holder inequality with 

/ = A\ (40) 
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g = A^, 

k — m + 1, 

k' = 1 + 1/m, 
drd0 

ax = 



sin^ ' 



we have 



where m 



(/ 

\Jr>l 



.drd9\ ™+i 



sm 



r>l 



r>l 



drdO 



(41) 
(42) 
(43) 

(44) 
(45) 



and is an integer greater than 1. With the power-law definition of Q in 



equation (18), it is easy to show that 



7 



4(m + 1) 



7 



drde 
sin 9 
drde 



H 



R 



An 



m + 1 jj.>i 



7 



sin 9 ' 



m + 



lr>l 



Now inequahty (45) leads to 



Hr < AnF^+' E;^+' 



16(m + 1^ 
7 



sin^ 



1 



2m+2 



(46) 
(47) 
(48) 



(49) 



Define 



H 



y = 



R 



16(m+l) 

7 



2m+2 



(50) 



Inequality (49) is then equivalent to say that y < 1. 

Figure 4 presents the y vahics, plotted against F^. for the families of n = 5, 7, 9 fields, 
using numerical solutions in Flyer et al. (2004) again. We see that the y values for these 
three families of fields vary between 0.75 to 0.95, with a tendency for the overall y values 
to increase with the n values. The fact that these y values are all less than 1 is just the 
verification of inequality (49), which is rigorous for all power- law axisymmetric force- free 
fields. The more interesting point is that these y values can get close to 1. This implies that 
we are getting close to finding a stringent bound on the total magnetic helicity by inequality 
(49), a more stringent inequality than inequality (35) is. 
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As pointed out by the referee, an interesting feature from Figure 4 is that our y values do 
not vanish to zero when both the total magnetic helicity and azimuthal flux vanish. Inserting 
equations (46)- (48) into equation (50), we find that our y values can be represented as 



y 



■I >» + 2 drdd 
> 1 bill d 



(51) 



which shows that values of y are not directly dependent on 7. Applying potential field 
^ = we get y = 0.88,0.90,0.92 for n = 5,7,9 respectively. These numbers are 

consistent with the tendencies of y values in Figure 4. 

Now multiply across equation (20) with flux function A and then integrate it with the 
element to get 

siny o 



f {m + 2)Q'dr^ = f 
Jr>i smfc/ J 



A— 

dr 



sin 6 



^jjBl^Bl^BDdr^^ 



(52) 

— - Jr>l ' ■ • Siiit7 

The left-hand side of the equation is proportional to E^p^ whereas, on the right hand side, 
the flrst term is bounded by the Chandrasekhar Virial Theorem and the second term is 
proportional to the magnetic energy Eq. Direct derivation then gives 



(m + 2)E^ < -Emax + Eo , 



(53) 



where Emax is given by equation (4). 



As Eq < Emax and Emax IS bouuded by the boundary condition at r = 1, the above 
inequality (53) shows that {171 + l)Ey, is also bounded. In consequence, we have E^ ^ 
as the free index m — > 00. This means that as we go to higher and higher n, the energy 
contributed by B'^ becomes neghgible. Much of the energy in the sum Bg + B'^ is in the 9 
component. In physical terms, we conclude that for the large-n families of fields, the fiux 
rope in the field is wounded so tightly that it approaches the highly localized structure of 
a line current or a sheet current outside of which B^^ is negligible. In other words, the free 
energy is stored as Bq in this exterior region associated with the line or sheet current. 



Applying inequality (53) to inequality (49), we get 



Hr < A7rF^+' 



16(m + l){^Emax 



En 



{m + 2)2 7 



2m.+2 



(54) 



a tight bound on Hr. 

To study this inequahty, let us define 
"16(m+l)" 



2m+2 



(55) 
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16(m+l)(|£;^„, + £;o) 



2m+2 



-'-'max ~ -'^0 \ 



Ci ( ^ " ) . (56) 



(m + 2)2 7 
Inequality (54) now can be written as 

m 

Hr<AttF^+'C2 . (57) 

The top panel of Figure 5 shows the C2 values, plotted against the azimuthal fluxes, for 
n = 5, 7, 9 families of fields. We see that as n increases from n = 5 to n = 9, the modulation 
of C2 values decreases, with the C2 values approaching 1. This indicates that as n — cxd, C2 
values tend to 1. This can be understood by the following analysis. C2 is a multiplication of 
Ci and ^ sEmax+Eo ^ 2m+2^ Sluce Ejnax, Eq and 7 are all bounded, the term ^ sEmax+Eo ^ 2m+2 
goes to 1 as m — CX3 except for 7 = 0. When 7 = 0, H^i = 0. This is the trivial potential 
field and is not an interesting field for our consideration of helicity upper bound. Owing to 
its mathematical form in equation (55), Ci approaches 1 as m ^ 00. This is indicated by 
the plot in the bottom panel of Figure 5, where Ci values are plotted against m and are 
approaching to 1 as m increases. So, on the limit of n — > 00 or, equivalently, m —>■ 00, 

Ci ^ 1, 
C2 ^ 1, 

F^+' ^ , (58) 

and hence we have 

Hr < AnF^ . (59) 

Now we see that in the limit of n — > 00, the bound on Hr reduces to a problem of 
finding an upper bound on F^. The force-free solutions of Flyer et al. (2004) suggest that 
such an upper bound exists (see also Figure 2 in this paper). Our tighter inequality (59) has 
thus provided a more stringent bound on Hr that had eluded from the lose inequality (35). 

Figure 6 shows the variations of Hr/{At:F^) against 7 for n = 5, 7 and 9 families. These 
plots show that the relationship Hr < AttF^ we derived forn — > 00 is also valid forn = 5, 7, 9 
families. So this relationship may be a general relationship that is vafid for all n values. If 
this is true, it tells us that the upper bound on the total magnetic helicity is a less stringent 
upper bound than that on the total azimuthal fiux. 



4. Conclusions and Discussions 



In this paper we proposed a hydromagnetic conjecture that there is an upper bound 
on the total magnetic helicity that a force-free field in an unbounded domain can contain. 
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Observations have shown that accumulations of hehcity with a preferred sign takes places in 
the two solar hemispheres. If our conjecture is valid, the relentless accumulation of magnetic 
helicity in the corona will lead to a CME-type eruption as a natural and unavoidable product 
of coronal evolution. This conjecture deserves further investigation to find a mathematical 
rigorous proof, or disproof. For the present, evidence for and implications of this conjecture 
can be found in the several families of power-law axisymmetric force-free fields governed by 
equation (20). 

We derived two rigorous inequalities (SI and S2 in Table 1) by applying the Cauchy- 
Schwartz and Holder inequalities, respectively. The former (SI in Table 1 or Equation (35) in 
the text) puts an absolute upper bound on the magnitude of total magnetic helicity for each 
fixed index n. This inequality is simple and rigorous but it is not tight for our purpose, with 
the bounds increasing monotonically with n. The second inequahty (S2 in Table 1) relates 
total magnetic helicity, azimuthal fiux and magnetic energy together in an interesting way, 
which presents a more stringent bound than SI does. 

Inequality S2 is further reduced to the simple form of Inequality S3 in Table 1 when 
n goes to infinity. We have not proved that the simple form of Inequality S3 is valid for 
all n although the numerical data of Figure 6 suggest it is. The vahdation of Inequahty 
S3 for large n is sufficient for our purpose, which is to remove our main concern arising 
from Inequality SI of whether the total helicity is still bounded when n becomes very large. 
Inequality S3 together with the numerical results from Flyer et al. (2004) suggest that such 
an upper bound on the total magnetic helicity exists even when n increases without bound. 

An important result of our study is the simple formula (SO in Table 1) for the relative 
helicity of an axisymmetric field. This formula is useful to future studies of axisymmetric 
fields applicable, independent of whether the field is force-free or in force balance. It has a 
simple physical interpretation: The total relative magnetic helicity in axisymmetric fields is 
just a simple convolution of the local wrapping of the azimuthal fiux element B^rdrdO with 
the poloidal flux A. 

We note that in Hu et al. (1997) the right hand side of Equation SO has been defined to 
be a form of helicity and the authors have showed that this form of helicity is conserved under 
the ideal hydromagnetic induction equation. Their derivation was carried out completely 
independent of the concept of relative helicity introduced by Berger & Field (1984). However, 
our derivation of Equation SO here, from a direct application of the concept of the Berger 
relative hehcity, shows that the helicity defined by Hu et al. (1997) is actually the same 
physical quantity as the Berger relative helicity for an axisymmetric field. 

What is the relationship between magnetic helicity accumulation and free magnetic 
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energy storage? Which one plays a more fundamental role in producing CMEs? It seems 
clear that the greater the total helicity the greater the magnetic energy would be, as suggested 
by the inequality (2). Although this inequality is rigorous only for a finite domain, a similar 
conclusion is expected generally in the unbounded domain that higher twist implies greater 
magnetic energy. The fact that the magnetic energy of a force-free field in the unbounded 
domain has an upper bound suggests that our conjecture on a similar kind of bound on the 
total helicity is probably physically sensible. It is interesting to observe from our Figure 2 
that the maximum storages of magnetic helicity and magnetic energy do not occur at the 
same values of 7. We suggest that magnetic helicity as a conserved quantity is for sure to 
be accumulated in the corona whereas the corona may release its free energy by magnetic 
rcconncction without leading to a CME expulsion. So for fields with large energy storage 
but moderate helicity storage, the field may release energy when a suitable trigger acts on it. 
But when magnetic helicity has accumulated to cross over the upper bound, we proposed, a 
CME expulsion becomes unavoidable. 

In this sense, exceeding the allowable level of helicity in a force-free field is only a 
sufficient condition for an eruption. It is an extreme condition where a field has accumulated 
enough helicity that an eruption becomes unavoidable. A CME expulsion may still occur 
even before this helicity limit is reached, as long as the separate necessary condition, that 
of having enough magnetic energy for an eruption, is met. For example, mass-loading is 
an efficient way to store enough free magnetic energy to drive CMEs, without having to 
cross the threshold on excessive magnetic helicity accumulation. In Zhang & Low (2004) 
an equilibrium state is given where prominences are supported by the magnetic fields in the 
so-called normal configuration. These solutions show that significant magnetic energy can 
be stored with only a relatively small amount of total azimuthal magnetic fiux and magnetic 
helicity. An interesting question for future investigation is then whether a coronal field may 
accumulate helicity to the point of exceeding the applicable upper bound without first having 
enough free magnetic energy to erupt into a CME. 
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Table 1. Significant equation and inequalities derived in the paper 



Number 



Equation or Inequality 



No. in the paper 



SO 

SI -Stt 

S2 
S3 



Hr = 47r X->i A X B^rdrde 



n+3 
2 



^0 + ^co' 



<Hii< 871 



n+3 
2 



(m+2)2 7 



2m+2 



(Equ. 17) 
(Equ. 35) 

(Equ. 54) 
(Equ. 59) 
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Fig. 1. — Shown is an axisymnictric magnetic field with its hncs of force seen projected on a r — 9 
plane, taken from Low (2001). The field lines are directed out of the shown plane if ^ 0. The 
middle panel shows an azimuthal field with Bj. = on r = 1, which cannot find an equilibrium. 
The upper panel shows a potential field with Br ^ on r = 1, meeting the necessary condition for 
an equilibrium. The bottom panel shows a field also meeting the necessary condition but with a 
detached flux rope, represented by the closed loops of poloidal projected lines. This configuration 
can be regarded as a sum of the above two configurations, of a field containing a certain amount 
of magnetic helicity and finding an equilibrium. 
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n=5 n=5 n=5 
0.601 I I 




Fig. 2. — Variation of the total magnetic energy (Eq), azimuthal flux (F^) and total magnetic 
helicity (Hji) vs. 7 along the solution curve for n = 5 (top panels), n = 7 (middle panels) and 
n = 9 (bottom panels). Each point, denoted by a plus symbol in the figure, represents a solution 
to equation (20). These curves of the solutions suggest that there may be an upper bound on the 
total magnetic helicity in addition to those on the total magnetic energy and azimuthal flux. See 
Flyer et al. (2004) for details of the solutions and the numerical solver. 
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n = 9 11 = 9 




Fig. 3. — Diagrams of the total magnetic helicity Hr (left panels) or of the total magnetic helicity 
divided by the total magnetic energy Hh/Eq (right panels) of each solution against the azimuthal 
flux F^p. Here Hr has been normalized by a factor of 47r. See text for discussions. 
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Fig. 4. — Plots of y values against azimuthal fluxes for n = 5, 7, 9 families of solutions. See 
text for definition and discussions. 
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1.30 F 




Fig. 5. — Top panel: Plots of C2 values against azimuthal fluxes for n = 5, 7, 9 families of solutions. 
See text for definitions and discussions. Bottom panel: Variation of Ci values against m. This plot 
shows that as m goes to infinity Ci goes to 1. 
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Fig. 6. — Plots of Hfi/ (AnFip) against 7 for n = 5, 7, 9 families of solutions. These plots show that 
the relationship Hr < AttF^ also valid for n = 5, 7, 9 families, so it may be a general relationship 
valid for all n values. 



